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Abstract 



We introduce a Hamiltonian for two interacting su(2) spins. This model is an analogue of 
the Bethe ansatz solvable p + ip-pairing Hamiltonian in that it admits a duality, from which the 
Read-Green and Moore-Read lines are present in the model's phase diagram. We use a mean- 
field analysis and exact Bethe ansatz results to investigate the ground-state properties of the 
system in the thermodynamic limit, denned here as the limit of infinite spin (or highest weight). 
We also numerically study the energy gap, fidelity and entanglement. We find that certain 
ground-state properties in the thermodynamic limit are ill-defined due to issues concerning 
£\j ' non-commuting limits of the system's coupling parameters and the thermodynamic limit. We 

^ , comment on possible implications for the p + ip-pairing Hamiltonian. 

. 1 Introduction 

rn 

\Q \ The p-Hp-pairing Hamiltonian has come to some prominence through the work of Read and Green pQ 

who, through a mean-field analysis, exposed a change in the topological properties of the ground- 
state wavefunction of the Hamiltonian as a function of the coupling parameter. An exact Bethe 
ansatz solution of the Hamiltonian was provided in [2], and shortly after two detailed analyses 
were conducted in [3113]. However their conclusions were different, with [3] arguing that there are 
three ground-state phases of the system while [3] maintained that there are two. In both cases 
the author's singled out the Moore-Read and Read-Green lines of the system, which will be defined 
later, as displaying some unique properties. In [3] it was argued that both lines should be identified 
with quantum phase boundaries, whereas [4] claimed that only the Read-Green line represented a 
quantum phase boundary. 

We consider a simplified Hamiltonian for two interacting su(2) spins, which we present as an 
analogue of the p + ip-pairing Hamiltonian in that it displays both the Moore-Read and Read- 
Green lines in the phase diagram as a result of a duality relation. We then conduct an analysis 
of the system using a mean-field approximation and the exact Bethe ansatz solution to investigate 
ground-state properties in the thermodynamic limit, defined here as the limit of infinite spin (or 
highest weight). An advantage of dealing with this Hamiltonian of two interacting su(2) spins is 
that analytic results are accessible. We also numerically solve the system to study the ground-state 
energy gap, fidelity [SJ[3], and entanglement .7,81, which have been used as indicators of quantum 
phase transitions in other simple bipartite quantum systems e.g. [SHU] ■ We find that certain ground- 
state properties in the thermodynamic limit are ill-defined, due to issues concerning non-commuting 
limits of the system's coupling parameters and the thermodynamic limit. We comment on some 
possible implications for the p + ip-pairing Hamiltonian. 
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2 The Hamiltonian 

The su(2) algebra has commutation relations 

[S+, S~] = 2S Z , [S z , S ,± ] = S ± . (1) 

Given a finite-dimensional module with highest weight s the action of the generators on the weight 
basis is 

S z \s,m) — m\s,m), 

S + \s, m) — y/ (s — m)(s + m + l)|s, m + 1), 
S^ls^m) = y/ (s + m)(s — m + l)|s, m — 1) 

with m = — s, — s + 1, s — 1, s. 

We will study the following coupled-spin Hamiltonian 

= (1 + G)(4{St + s) + 4(S% + a)) - G{ Zl S+ + z 2 S+){z 1 S{ + z 2 S^) (2) 

which acts on the (2s + l) 2 -dimensional tensor product space with basis 

{|s, mi) (g) |s, m 2 ) : mj — —s, —s + 1, s — 1, s; j = 1, 2}. (3) 

Throughout we assume zi > z 2 > 0. For given s the model has three coupling parameters. One 
of these can be fixed by choosing the overall energy scale, e.g. setting z\ = 1, which leaves two 
independent parameters z 2 and G. Define 



H 


= 4W + 


si) + zl{S z + s 2 


Q + 


= ziSt + 


z 2 S£ 


Q~ 


— + 


z 2 S 2 


B+ 


= -S+ + 






Zl 


z 2 2 


S z 


= si + si 





such that 

H = (1 + G)H - GQ+Q- (4) 
and [S z , H] = 0. From the commutation relations (JT]) the following equation can be shown to hold 

[H, [B+) M ] = MQ+(B+) M - 1 (1 + GM + 2GS Z ). 

If a state has spin m (the eigenvalue of S z ) the above means that choosing 1 + GM + 2Gm = 0, so 
M = —2m — 1/G, the two states and = (B + ) M have the same energy. We call the states 
\ip) and dual states. Since the state has spin m the state l^) has spin m + M = — m — 1/G. 
We identify the following cases: 

• The dual of the lowest weight state = |s, — s) ® |s, — s) occurs when M = 4s — 1/G leading 
to 1^) having spin m = 2s — 1/G. This is referred to as the Moore- Read line, on which the 
ground-state energy is zero. 

• If M = then a state is self-dual when m = — 1/(2G). This is referred to as the Read-Green 
line. 

Note that the above definitions for the Moore-Read and Read-Green lines are independent of the 
parameters Z\ and z 2 , analogous to the case of the p + ip-pairing Hamiltonian where these lines 
are independent of the single particle momentum distribution [2H4]. The Hamiltonian ([2]) may be 
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viewed as a degenerate, two-level limit of the p + ip-pairing Hamiltonian on a restricted Hilbert 
space, which has been studied in [12] . 

Our main objective is to investigate ground-state properties of the model in the thermodynamic 
limit which we define by a sequence of ordered triples (G,m,s) such that G — > 0, m — ^oo,s^oo 
with p = m/(2s), g = AsG held constant. In terms of these parameters the Moore- Read line is given 
by p = \ — 2/g while the Read-Green line is p = —1/g. 



3 Mean-field analysis 

The initial study by Read and Green [T] on the p + ip-pairing Hamiltonian was undertaken using a 
mean-field analysis. Here we start by taking an analogous study for the coupled-spin model. The 
Hamiltonian can be approximated by a Hamiltonian with linear terms through 

AB « A(B) + (A)B - (A)(B) 

where (A) denotes the expectation value of the operator A. The mean-field approximation which 
linearises the Hamiltonian ([2|) is 

H m{ = (1 + G)s(zf + z\) + pm + ((1 + G)z\ - p)S{ + ((1 + G)z 2 2 - p)S* 
A A* A 2 

where p is a Lagrange multiplier (the analogue of the chemical potential) which is introduced as the 
mean- field approximation breaks the conservation of S z , and 

A = 2G(ziSr + z 2 S 2 -). 

Setting A = |A| the ground state energy is found to be 

E = (1 + G)s{zl + 4) + pm +^- S ^{{l + G)zl-pf + zi& - s^j ((1 + G)z\ - p) 2 + zfA 2 . 

Using the Hellmann-Feynman theorem we can determine the mean- field expectation values: 

/ dH mi \ = £)Eo 
\ dG I dG 



A^_ = a (l + G)z 2 -p + A 2 /G ^ (1 + G)z 2 - P + A 2 /G 
2G 2 SZl ^((l + G)z 2 -p) 2 + z 2 A 2 2 V((l + G)z 2 -p) 2 +z 2 A 2 ' 



Taking the thermodynamic limit gives 
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9 y/(z( - pf + z 2 A 2 ^(z$ - p) 2 + z 2 A 2 



(5) 



Also 



dH mf \ dEn 



dp I dp 

s((l + G)z\ - p) + G)z\ - p) 



-(SI) - (Sf) = 



v/((l + G)z{ - pf + z 2 A 2 + G)zi - p) 2 + z|A 2 

Taking the thermodynamic limit and using ([5]) gives 

g + 2 M = - P P (6) 

9 ^{z{-pY+zi^ y/(zi - pY + z 2 2 A 2 
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Returning to the energy expression in the limit s — > oo, and taking into account Eqs. (|5l6p . we 
obtain the intensive energy expression 

,. Eq 
eo = hm — 



= (z( + 4) + 2pp+ — -^(z 2 - p) 2 + z 2 A 2 y/(4 p) 2 + z 2 A 2 
= 2 / _ 1 2z 2 + A 2 -2p \ J 1 2z 2 + A 2 -2p \ 

X \ 2 y / {z f-p)2 + z - tA 2) 2 [ 2 ^_ F _ f _y. 

It may be verified that setting p = A 2 /4 in (|5I6[) and combining these equations leads to the equation 
of the Moore- Read line. Furthermore, setting p = A 2 /4 in ([7]) gives eo = in agreement with the 
result that ground states on the Moore-Read line are dual to the lowest weight state. 
Recalling the assumption z\ > Z2 the mimimum excitation energy (i.e. the gap) is 



\j{zl-p) 2 



AE = ^/ '(z 2 - p) 2 + z 2 A 2 



which vanishes when z 2 = p = 0. However turning to ^ we see that the limits z 2 — > and p — s- 
do not commute in that equation. First setting p = leads to the condition of the Read- Green line. 
Alternatively first setting z 2 = yields 

1 p 1 , . 

V = o llm n ] \ 8 ) 

2 p^O \p\ g 

which shows a discontinuity as p passes through zero. This indicates that while gapless excitations 
occur in the limit z 2 ,p —> 0, where they are to be found in the phase diagram will depend on how 
these two limits are taken (at least within the mean- field approximation). 

In the analogous study of Read and Green [T] , a similar type of discontinuity was also observed 
in a particular number operator expectation value. Through this observation the Read-Green line 
was identified with a topological phase transition. For the simplified coupled-spin model we will 
next conduct an analysis using exact results. In this manner we will highlight some subtle aspects 
regarding the ground-state properties of the system, which appear to be beyond the scope of a 
mean-field analysis. 

4 Exact results 

4.1 The limiting case z 2 = 

Here we set z 2 = from the outset, and fix the energy scale by letting z\ — 1, in which case the 
Hamiltonian is simply 

H = (1 + G)(ST + s) - GS+S^. 

The above Hamiltonian is diagonal in the basis ([3]) with eigenvalues 

E = mi + s - G(s 2 - m\). (9) 

For each value of m — mi + m 2 we can determine the ground state energy, and the gap, by 
considering where level crossings occur. We need to consider several different cases. Starting with 
m = mi + m 2 < we have 

mi = — s, m 2 = m + s =>■ E = 0, 

mi = -s + l, m 2 = m + s-l =s> E = (1 + G) - 2sG. 
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Figure 1: Phase diagram obtained from duality. The solid line fj, = 1 — 2/g is the Moore-Read line, 
on which the intensive ground-state energy is zero. The dotted line \i = — 1/g is the Read-Green 
line. The mean-field analysis predicts that there are gapless excitations on the Read-Green line as 
Z2 — > 0, but also that the limits /i — > —1/g and Z2 — > do not commute. With respect to the 
duality discussions of Sect. 2, the Moore-Read line denotes states dual to the lowest weight state 
(/j = —1), while the Read-Green line denotes states which are self-dual. Note that neither boundary 
line depends on the variables z\ and z 2 - 



In the thermodynamic limit it is seen that that for to < a first level crossing occurs at 

3 = 2. (10) 

Next consider 

m 1 =s + m-l, m 2 = -s + l => E = 2s + to - 1 + G((m - l) 2 + 2s(m - 1)), 
mi = s + m, m2 = — s ==>■ E = 2s + to + G(m 2 + 2sm). 

Here there is an energy gap 

AE = |l + G(2m- l + 2s)|. 
In the thermodynamic limit the final level crossing occurs when 



which requires that /i < —1/2. 

Now we consider the case for m > 0: 



mi = — s + to, to 2 = s =>■ E = to + G(m 2 — 2sm), 

toi = -s + to + 1, m 2 = s-l =^> E = to + 1 + G((m + l) 2 - 2s(m + 1)). 

The gap 

A£;= |1 + G(2to+1-2s)| 
goes to zero in the thermodynamic limit when 

"=H- <12) 
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This requires that < fi < 1/2. In the region bounded by the equations (|10I11I12[) the level crossings 
become dense in the thermodynamic limit, giving rise to a gapless phase. The following relations 
hold for the gap between the ground and first-excited states in the regions depicted in the phase 
diagram of Fig. 2: 

(13) 
(14) 
(15) 
(16) 



I : 


AE = 


1 + 


9{fJ> ~ 


1/2), 


II : 


AE = 


1 - 


5/2, 




III : 


AE = 


-1 


-s(m 


+ 1/2 


IV : 


AE = 


0. 







1 




Figure 2: Ground-state phase diagram for z 2 — 0. Expressions for the gap between the ground and 
first-excited states in the thermodynamic limit are given by Eqs. (|13H16p . Region IV is gapless. 



There is clear discrepancy between Figs. 1 and 2, in particular for the predicted gapless regions. 
In deriving the two phase diagrams we consider behaviours ass— > oo, m — >■ oo, G — > with 
g = 4sG, fi = m/(2s), and as z 2 — > 0, yielding phase boundary lines as functions of the coupling 
parameter g. So far there are indications that the limiting behaviours do not always commute. 
Similar issues with non-commuting limits have been observed in the different context of attractive 
bosons [13]. In the following three Subsections we will appeal to further exact results to gain a 
better insight into the properties of the system. 



4.2 The case m = 1 — 2s 

In the sector where m = 1 — 2s the Hilbert space is two-dimensional and exact results are easily 
extracted. In matrix form, with respect to the basis 

{\s, -s + l)<8) \s, -s), \s, -s) «) \s, -s + 1)}, 

the Hamiltonian is 

tt _ ( (1 + G - 2sG)z\ -2sGz lZ2 \ n? , 
^ -2sGz lZ2 (l + G-2sG)zf J ( > 

with the ground-state energy 

E [m ^ s] - (1+g/(45) 7 /2)(z? + z|) - y{i+ 9 i&)- g im*i-w+9**& (18) 



G 



and energy gap to the excited state 

A£ = ^(1 + g/(M - g/2) 2 (z 2 - z 2 ) 2 + g 2 z\z 2 . 

Here fi — > — 1 in the thermodynamic limit, which at g = 1 corresponds to the limit where the 
Moore-Read and Read-Green lines meet. It is found that 

lim lim lim AE = 

22—^0 g— >2 s— >oo 

with the result being invariant with respect to interchanging the order of the limits. Moreover, as 
s — > oo the only solution for AE = when Z\ ^ is given by z 2 = 0, g — 2, meaning that gapless 
excitations do not occur at the Read-Green line in this limit /x — > — 1. The behaviour accords with 
Fig. 2 rather than Fig. 1. 
Next consider 



8{AE) g{zl + z 2 ) 2 - (2 - (1 - g + g / '(4s)) / ' s)(z 2 - z 2 2 ) 2 



We find 



<9# 



lim lim lim 

22—^0 g— >2 s— t-oo 

lim lim lim 

</^2 + 2 2 ^0 S^OC 

lim lim lim 

g^2~ 2 2 — >0 s— >oo 

lim lim lim 

s— >oo 22— >0 g— >2 



4A£ 



a(A^) 


= 0, 


dg 


d(AE) 


7 2 
Z l 


dg 


~ 2 ' 


d(AE) 


z 2 
z i 


dg 


2 


d(AE) 


Z l 


dg 


2 



lim lim lim — ^ — - 

z 2 -t0 s^mj->2 

lim lim lim — ^— — - 

3^2+ s^oo z 2 ->0 

lim lim lim — ^— — - 

g^2- s^oo 22-i-O 

lim lim lim — - 



z l 



This establishes that the non-commutativity of limits is not just an artefact of the mean-field ap- 
proximation, but is an intrinsic property of the model. 

4.3 Non-analytic behaviour on the Moore-Read line 

Next we will use the duality properties of the model to show that non-analytic behaviour is found 
on the Moore-Read line appearing in Fig. 1. Our approach is through a study of how the intensive 
ground-state energy e (g) varies with g while zi, z 2 and fi are held constant. 

First we consider two ground states of the system at different values of the coupling parameter: 



G« - 



1 



2s — m 



G< 2 > = 



1 



2s — m — 1 



such that AG = G( 2 ) - G« = G^G<- 2 \ We also have Ag = 4s(G< 2 ) - G«) = 4sG«G< 2 ). It is 
found that in the thermodynamic limit 



lim = lim 4sG (1) 

s— >oo s— >oo 

lim g {2) = lim 4sG (2) 



lim 

8— >00 

lim 



4s 



s^oo 2s — m 
4s 



= lim 



so we hereafter set 



Then 



>°o 2s — m — 1 
2 

3mr = z • 

1 -/i 



lim 



2 



oo 1 -n- l/(2s) 1-/J, 



lim 4s Ag = 5mr . 

s— too 



(19) 
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The ground state corresponding to = 4sGW is on the Moore-Read line, and so has energy 
^oGr 1 ') = 0. The ground state corresponding to — AsG^ is dual to the ground state of (fT7| . 
and so has energy E {g^) — Ey m= i_ 2s \ as given by (fT5|) . We also have 

eo(5MRj = hm = hm = 0. 



Now consider in general terms 



de (ff) _ e (ff + 6)- e (g) 
dg <5->o 5 

= lim lim ML+iL^M. 

5— >0 s— >oo so 



One could alternatively define 

de (g) = Hm l dE (g) 
dg s-s-oo s dg 

,. ,. ^(g + Agj-i^g) 
= lim lim . 

s— s-oc <5— >0 so 

The definitions will be equivalent if the limits S — > and s — ¥ oo commute. 
If the limits g — > gyiR and s — > oo commute we have 

deo(ff M R.) _ e (g M R + S) - e (g M R) 
dg <5-s-o S 

- V e o(ffMR + $) 
<5->0 5 

E (gW + 5) 



lim lim 
lim 



<5->0 s— >oo sS 

£ (g (1) +Ag) 



sAg 

lim — lim S (g (2) ) 

(2 = 9M f Zl + - - J-y/fi W2) 2 (*? - z\Y + g^zl. 



9mr 9mr 
where in the last step we use (|18I22[) . Alternatively 



de (g M R) _ j. e (g M R + 6) - e (g M R) 
dg s^o S 

- lim e °( ffMR + ^ 



s^o S 
E (g^ + S) 

5—>0 s— >oo s5 

Eo(gW-Ag) 



lim lim 

<S->0 s- 

— lim 



sAg 



= - lim -— lim E (g {1) ) 

s— too sA(J s— too 

= 0. 

So assuming that the limits s — > oo and g — > gMR commute leads to the partial derivative being 
discontimuous. 



8 



4.4 Bethe ansatz solution 



As in the case of the p + ip model [5] , the Hamiltonian admits a Bethe ansatz solution. The 
energy spectrum is given by 

M 

E = (i + G)Y,yk 

fc=i 

where the y^, k = 1, M, are roots of the Bethe ansatz equations 

M 

^ E A (20) 



■ 2M - 4s - 1 



2s 



2s 



;(/*■ 



2/A 



2/A 



, Vk - Vj 



and the total spin is m = M — 2s. For each solution of the Bethe ansatz equations the associated 
eigenstate is 

M 

\i>) = l[C(y k )(\s,-s)®\s,-s)) 



where 



fc=i 



C(y) = 



ziS+ 



y- 4 



(21) 



The validity of this solution can be verfied in a direct manner, using only the commutation relations 
(HJ). We omit the details. 

Following [3,4 , we look to solve the Bethe ansatz equations in the thermodynamic limit under 
the assumption that in this limit the roots of (|2"D)l are densely distributed on a curve V in the 
complex plane which is invariant under reflections about the real axis. Introducing r(y) as the 
density function on T for the roots, we divide (|20[) by 2s and take the limit s —¥ oo to obtain 



2g~ x + 2/i 



1 



y - z\ y- z% 



= 2P / \dy 



, r(y>) 
y-y' 



(22) 



y y~*i 

where P refers to the Cauchy principal value for the integral. To solve the problem we define an 
analytic function h(y') outside T such that 



r(y') = h + (y')-h-(y') 



(23) 



where h+(y') and h-(y') are the limiting values of h{y') to the right and left of T. We can then 
rewrite (|2"21 as 



2g- 1 +2/i 1 



1 



y 



y-4 



y 



= 2 <b dy 



j Hy>) 
y-y' 



where L encircles T. Next we take the following ansatz: 



R(y) 

2m 



2i 



^2 



S 



y y 



(24) 



(25) 



where 



R{y) = y/(y-a)(y-b) 

with a, b = a* denoting the end points of the arc T. Replacing (1251) in 
integral by residues leads to the conditions 



(26) 

and evaluting the contour 



S 

s 



2 v /( Z f-a)(z2_6) : 
<j>{zi) + 0(z 2 ), 



n b 
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These combine to give 



l+2»= - Va b - + - j ± (27) 

9 * y/(zi-a)(zf-b) y/(zi-a)(zi-b) 



We may also compute 



/* = -l + J W\r{y') 
= -1+ I dy'h(y') 



L 



a - 



S - - z 2 <7V2) s S 



9 y/(z*-a)(z*-b) y/(z*-a)(z*-b) 



(28) 



Identifying p = ^/ab and A 2 = a + b + 2^/ab, Eqs. ()27|28|) are precisely the mean-field equa- 
tions ([6 1 "ill respectively. This would suggest that the mean-field equations should be exact in the 
thermodynamic limit as they are reproducible from the exact Bethe ansatz solution. However, the 
above calculations assume that in the thermodynamic limit the roots of the Bethe ansatz equations 
become dense on a curve T. This assumption is not always valid. Suppose that the curve T closes 
such that the end points a, b become real and equal. Then Eqs. (|27l28p reduce to 



2 

9 ' _f ~ \ z l ~ a \ ' \ Z 2 



^=TZT^ + TJ^y (29) 
2 _ n \ ■ 



9 \ z l - a \ k 



(30) 
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If a > the only solution is /1 = — 1, whereas if a < we obtain the equation of the Moore-Read 
line 11 = 1 — 2<7 _1 . But the ground-state roots on the Moore- Read line are all zero, as the ground 
state is dual to the lowest weight state. Thus taking s — > 00 after g — > <?mr leads to a delta 
function distribution of roots at the origin, while taking g — > g^j R after s — > 00 leads to the roots 
being distributed on a closed curve enclosing the origin. The non-commutativity of these limits has 
previously been observed in the p + ip-pairing Hamiltonian. In 2,3 this was referred to as a zeroth 
order transition. We may identify the Moore-Read line as the boundary at which mean-field theory 
breaks down. In the region immediately beyond the Moore-Read line we cannot assume that the 
ground-state roots of the Bethe ansatz equations are densely distributed on a single curve, in order 
to produce the mean-field equations from the Bethe ansatz solution. 



5 Indicators of quantum phase transitions 

While there have been numerous studies of quantum phase transitions, there still lacks a definitive 
property which may be used to identify them. In the remaining Subsections we study various 
methods to gain an insight into the nature of possible quantum phase transitions for this system. 
Using numerical diagonalisation of the Hamiltonian ([2]), below we compute the energy gap, fidelity, 
and entanglement (as given by the von Neumann entropy). Each quantity is plotted as a function 
of the parameter g~ x . We fix fi = —0.99, Z\ — 1, and adopt the parameterisation z-i — (2000s) -1 / 2 
such that Z2 as s — > 00. We focus on this case as it permits us to compactly cover the regions 
II, III and IV in Fig. 2. 
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5.1 Energy Gap 



In Fig. [3] we plot the energy gap. There are two different points of minima corresponding to the 
boundary lines between region III and IV, and between region IV and II, in Fig. [2] It is seen in the 
inset that as s increases the energy gap decreases in region IV of Fig. 2. The results appear more 
consistent with Fig. 2 rather than Fig. 1, as to be expected in view of the results from Subsect. 4.2. 




Figure 3: Energy gap between the ground state and the first excited state as a function of g , with 
/i = —0.99. The results shown in the inset indicate the tendency towards vanishing gap in region 
IV as s — > oo. The solid black line represents the analytic curve for the energy gap as given by Eqs. 
flUHS} for when z 2 = 0. 



5.2 Fidelity 

Another means to investigate quantum phase transitions is using fidelity 5,6 , a concept widely used 
in quantum information theory. The fidelity of two states is simply defined as the modulus of the 
inner product 

7-ftM) = |W#>l- (31) 

Computing this quantity for ground-state vectors associated with Hamiltonians differing only by 
a small change in the coupling parameter, the fidelity may be sensitive to variations of the cou- 
pling parameter. Generally a minimum in the fidelity points to a rapid change in the ground-state 
properties, which may indicate a critical point of the system. 
Here we take two coupling values 

{g ± )- 1 =g-\l± 1 ) 

with 7 = 0.001 to compute the fidelity of the ground states associated with g ± 1 and then consider 
the behaviour of the fidelity as g is varied. As is seen in Fig. IH for values of the parameter g^ 1 
above 0.5 and values below 0.49, the fidelity is close to 1, but suffers an abrupt change at these two 
points. The two points correspond to the boundary lines between regions II and IV and regions IV 
and III for [i — —0.99 in the phase diagram Fig. [5J Also, it can be seen that the minimum of the 
fidelity is decreasing as s increases. 
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0.54 



Figure 4: Ground-state fidelity as a function of g^ 1 , with /i = —0.99 and 7 = 0.001 for different 
values of s as indicated in the legend. The results show significantly decreasing fidelity in region IV 
as s — > 00. 



5.3 Entanglement 

The coupled-spin model is naturally seen as a bipartite system. In this case we can use the von 
Neumann entropy of the reduced density operator to quantify the entanglement of the system [7J[9]- 
[TTj . Using the fact that the total spin m is conserved, we can write a general state of the system 
(for m < 0) in terms of ((3|) by 

m+s 

1*)= J2 d j \s,j)®\s,m-j), (32) 

3=—a 

for some complex numbers dj . The density operator for state (|32l) is given by 

m+s 

p=|*)(*|= d i d* j \s ) i)(s > j\^\s,m-i)(s,m-j\. (33) 

i,j=-s 

We take the partial trace with respect to the spin sector 2 yields the reduced density operator for 
the spin sector 1: 

m+s 

Pi=tr 2 (p)= J2 \dj\ 2 \s,j)(s,j\. (34) 

j=—s 

The von Neumann entropy of entanglement for the state is given by 

m+s 

£{ Pl ) = -tr[ Pl log 2 ( Pl )] = -J2 l^| 2 log 2 (M,| 2 ). (35) 

j=—» 

Using the formula (|35p we compute the entanglement of the ground state with fi = —0.99 as function 
of g^ 1 , as shown in Fig. [5j We can identify two abrupt changes in the entanglement value about 
the points g^ 1 — 0.49 and g _1 = 0.5 in agreement with energy gap and fidelity analysis. Note that 
the entanglement variation gets sharper as the value of s is increased. 
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Figure 5: Ground-state entanglement as function of g^ 1 for /z = —0.99 for different values of s 
as indicated in the legend. The results show significantly increasing entanglement in region IV as 
s — > oo. 



6 Conclusion 

The mean-field analysis of the model predicted the ground-state phase diagram shown in Fig. 1 
displaying the Moore-Read line, where the ground-state energy is zero, and Read-Green line, on 
which the gap goes to zero as z-i —> 0. These lines also appear naturally with regards to the duality 
described in Sect. 2 and are independent of the variables z\ and z%. This is exactly the same phase 
diagram as that found for the p + ip-pairing Hamiltonian [2]-[4]. However starting with a finite- 
size system, setting Z2 = 0, and then taking the thermodynamic limit gave the ground-state phase 
diagram shown in Fig. 2. In Sect. 4 we turned to the exact Bcthe ansatz solution for the model 
and through it found features associated with both phase diagrams. In the limit [i — > — 1 the gap 
was minimised at g — 2 as Z2 — > 0, in agreement with Fig. 2. Upon further analysis we also saw 
that the limits s — > oo and Z2 — > do not always commute, with the derivative of the energy gap 
being ill-defined in this limit as g — > 2. On the other hand the Moore-Read line of Fig. 1 was 
found to exhibit non-analytic behaviour if the limits s — > oo and g — ¥ <?mr were assumed to be 
commutative. Using the exact Bethe ansatz solution we argued that the Moore-Read line signifies 
the breakdown of the validity of mean-field solution. Finally we numerically studied the gap, fidelity, 
and entanglement to illustrate that signatures of the phase diagram in Fig. 2 were found in these 
quantities. 

The complexity observed here in the simplified couplcd-spin model highlights the need for deeper 
investigations into characterising the phase diagram of the p + ip-pairing Hamiltonian, in a manner 
which fully takes into account the peculiarities presented by the non-commuting limits. We expect 
that this is particularly the case in the so-called weak-pairing phase, which lies between the Moore- 
Read and Read-Green lines. It is in this region that the ground-state roots of the Bethe ansatz 
equations cannot be uniquely described by a density function, with support on a connected curve T, 
which varies smoothly with g. In this region there are limiting root distributions which have delta 
function terms at the origin, as a result of the duality, that are not connected to the distribution of 
the remaining roots. 
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